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Abstract. This paper presents a more complete version than hitherto published 
of our explanation of a transition from regular to irregular motions and more 
generally of the nature of a certain kind of deterministic chaos. To this end we 
introduced a simple model analogous to a three-body problem in the plane, whose 
general solution is obtained via quadratures all performed in terms of elementary 
functions. For some values of the coupling constants the system is isochronous 
and explicit formulas for the period of the solutions can be given. For other values, 
the motions are confined but feature aperiodic (in some sense chaotic) motions. 

This rich phenomenology can be understood in remarkable, quantitative detail 
in terms of travel on a certain (circular) path on the Riemann surfaces defined by 
the solutions of a related model considered as functions of a complex time. This 
model is meant to provide a paradigmatic first step towards a somewhat novel 
understanding of a certain kind of chaotic phenomena. 
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1. Introduction 

The fact that the distinction among integrable or nonintegrable behaviors of a 
dynamical system is somehow connected with the analytic structure of the solutions 
of the model under consideration as functions of the independent variable "time" 
(considered as a complex variable) is by no means a novel notion. It goes back to 
classical work by Carl Jacobi, Henri Poincare, Sophia Kowalevskaya, Paul Painleve 
and others. In recent times some of us had the good fortune to hear in several occasions 
such ideas clearly described by Martin Kruskal [221 123] ■ A simple-minded rendition of 
his teachings can be described as follows: for an evolution to be integrable, it should 
be expressible, at least in principle, via formulas that are not excessively multivalued 
in terms of the dependent variable, entailing that, to the extent this evolution is 
expressible by analytic functions of the dependent variable (considered as a complex 
variable), it might possess branch points, but it should not feature an infinity of them 
that is dense in the complex plane of the independent variable. 

A number of techniques collectively known as Painleve analysis have been 
resurrected and further developed over the last few decades (for a review see, for 
instance, [28j). In essence, they consider an ansatz of the local behaviour of 
a solution near a singularity in terms of a Laurent series, introducing it in the 
equations and determining the leading orders and resonances (terms in the expansion 
at which arbitrary constants appear). Painleve analysis has been extended to test 
for the presence of algebraic branching (weak Painleve property [2^ ) by considering a 
Puiseaux series instead of a Laurent series. These analytic techniques (which have been 
algorithmized and are now available in computer packages) constitute a useful tool in 
the investigation of integrability: in many nonlinear systems where no solution in 
closed form is known, Painleve analysis provides information on the type of branching 
featured by the general solution or by special classes of solutions. It has also proved 
useful to identify special values of the parameters for which generally chaotic systems 
such as Henon-Heiles or Lorenz are integrable [4l I16j . 

On the opposite side of the spectrum lie chaotic dynamical systems and it is 
natural to investigate the singularity structure of their solutions. Tabor and his 
collaborators initiated this study in the early eighties for the Lorenz system [211 and 
the Henon-Heiles Hamiltonian[16 . They realized that the singularities of the solutions 
in complex-time are important for the real-time evolution of the system. The complex 
time analytic structure was studied by extensions of the Painleve analysis involving 
the introduction of logarithmic terms in the expansion — the so called ^E'-series — 
which provides a local representation of the solutions in the neighbourhood of a 
singularity in the chaotic regime. Their local analytic approach was complemented by 
numerical techniques developed for finding the location of the singularities in complex 
time and determining the order of branching |14j . In all the chaotic systems under 
study, they observed numerically that the singularities in complex time cluster on a 
natural boundary with self-similar structure [15j . An analytic argument to explain 
the mechanism that leads to recursive singularity clustering was given in [24] . Similar 
studies relating singularity structure, chaos and integrability have been performed by 
Bountis and his collaborators. Going beyond the local techniques described above, the 
emphasis is put on a global property of the solutions: whether their Riemann surface 
has a finite or an infinite number of sheets. Bountis proposes to use the term integrable 
for the first case and non-integrable for the second [21 13] • Using mostly numerical 
evidence he conjectures that in the non-intcgrable cases the Riemann surfaces are 



Towards a Theory of Chaos Explained as Travel on Riemann Surfaces 



3 



infinitely-sheeted and the projection on the complex plane of the singularities is dense. 
Combining analytical and numerical results for a simple ODE, Bountis and Fokas 
[T5] have identified chaotic systems with the property that the singularities of their 
solutions are dense. 

Painleve analysis and its extensions are useful and widely applicable. However, as 
local techniques, they provide no information on the global properties of the Riemann 
surfaces of the solutions, such as: the number and location of the movable branch 
points a solution has, and moreover how the different sheets of the Riemann surface 
are connected together at those branch points. Understanding these global properties 
is important for the dynamics; a detailed analysis of the Riemann surface associated 
to the solutions of a dynamical system, whenever it can be done, provides a much 
deeper understanding than can be obtained by local techniques alone. 

This is precisely the motivation of the investigation reported herein: to introduce 
and study a model which is simple enough that a full description of its Riemann 
surface can be performed, yet complicated enough to feature a rich behaviour, possibly 
including irregular or chaotic characteristics. 

Such a model was initially presented in [1] and in this paper we continue 
investigating its properties. This line of research originates from a "trick" that 
is convenient to identify isochronous systems [7l [9] - a change of dependent and 
independent variables, with the new independent variable traveling on a path in the 
complex plane. Later it was shown that many isochronous systems can be written by 
a suitable modification of a large class of complex ODEs [TOl . Using local analysis 
and numerical integration in two many-body systems in the plane \lT\ . it was 
discovered that outside the isochrony region there exist periodic solutions with much 
higher periods as well as possibly aperiodic solutions, and the connection among this 
phenomenology and the analytic structure of the corresponding solutions as functions 
of complex time was illuminated. However, those systems were too complicated for a 
complete description of the Riemann surface to be achieved. 

Recent work along these lines includes problems whose solution is obtained by 
inversion of a hyperelliptic integral: the corresponding Riemann surfaces have been 
studied in [T71 \W\ , together with the implications on the dynamical properties of the 
models. 

In the present paper we provide many details that were reported without proof 
in [1], such as the description of the general solution by quadratures, and we also 
exhibit other properties of the model that were not present in [l], such as similarity 
solutions, equilibrium configurations and small oscillations. Our investigation of this 
model will continue in a subsequent publication, [12] , where the full description of the 
geometrical properties of the Riemann surface will be given. 

This paper is organized as follows: in Section 2 we present our model, including 
in particular the relationship among its physical version (independent variable: the 
real time t) and its auxiliary version (independent variable: a complex variable r), 
and we outline the main findings reported in this paper. In Section [3] we discuss 
the equilibrium configurations of our physical model, and the behavior of this system 
in the neighborhood of these solutions, and we also obtain certain exact similarity 
solutions of our model and discuss their stability. In Section |4] we discuss the analytic 
structure of the solutions of the auxiliary model via local analyses a la Painleve, since 
the analytic structure of these solutions plays a crucial role in determining the time 
evolution of our physical model. In Section [5] we show how the general solution of 
our model can be achieved by quadratures and in Section [6] we outline the behavior of 
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our model based on these results. Finally, in Section [7] we summarize our results and 
comment on future developments. This paper also contains a few Appendices, where 
certain calculations are confined (to avoid interrupting inconveniently the flow of the 
presentation) as well as certain additional findings. 

2. Presentation of the model 

In this section we introduce the model treated in this paper, and we outline our main 
findings that are then proven and further discussed in subsequent sections. 

2.1. The auxiliary model 

The auxiliary model on which we focus in this paper is characterized by the following 
system of three coupled nonlinear ODEs: 

^/ _ 9n+2 _|_ ffn+l 

Cn — Cn+1 Cn — Cn+2 

Notation: here and hereafter indices such as n, m range from 1 to 3 and are defined 
mod(3); r is the {complex) independent variable of this auxiliary model; the three 
functions Cn = Cn{T) are the dependent variables of this auxiliary model, and we 
assume them to be as well complex; an appended prime always denotes differentiation 
with respect to the argument of the functions it is appended to (here, of course, with 
respect to the complex variable r) ; and the three quantities g„ are arbitrary "coupling 
constants" (possibly also complex; but in this paper we restrict consideration mainly 
to the case with real coupling constants; this is in particular hereafter assumed in 
this section). In the following we will often focus on the "semisymmetrical case" 
characterized by the equality of two of the three coupling constants, say 

gi= 92= .9, 53 = /, (2) 

since in this case the treatment is simpler yet still adequate to exhibit most aspects 
of the phenomenology we are interested in. More special cases are the "fully 
symmetrical" , or "integrable" , one characterized by the equality of all three coupling 
constants, 

9 = f, 9i= 92= 93 = 9, (3) 
and the "two-body" one, with only one nonvanishing coupling constant, say 

5i = 32 = 5 = 0, 93 = f ^0 ■ (4a) 
In this latter case clearly 

C^=0, C3(r)=C3(0) (46) 
(see ([T])) and the remaining two-body problem is trivially solvable, 

Cs{r) = l [Ci(0) + C2(0)]-(-)^ |i [Ci(0)-C2(0)f + /r|'^', ,s = l,2, (4c) 

while the justification for labeling the fully symmetrical case ^ as "integrable" will be 
clear from the following (or see Section 2.3.4.1 of |5). Before introducing our physical 
model, let us note that the auxiliary system ([1]) is invariant under translations of both 
the independent variable r (indeed, it is autonomous) and the dependent variables 
Cn (t) {Cn (t) Cn (■'') + Coi Co = 0)j is morcovcr invariant under an appropriate 

simultaneous rescaling of the independent and the dependent variables. 
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2.2. The trick and the physical model 

The trick mentioned above, relating the auxiliary model to the physical model, 
amounts in our present case to the introduction of the (real) independent variable 
t ("physical time"), as well as the three (complex) dependent variables Zn = Zn{t), via 
the following positions: 

exp(2iujt) — 1 

T= TT- ' (5a) 

ZlUJ 

z„(i) = exp(-iwi) C«(t) . (56) 

We hereafter assume the constant uj to be real (for definiteness, positive, to > 0; note 
that ioT uj — the change of variables disappears), and we associate to it the period 

T=-. (5c) 

Note that this change of variables entails that the initial values z„ (0) of the "particle 
coordinates" z„(t) coincide with the initial values Cn(0) of the dependent variables of 
the auxiliary model ([1]): 

ZniO) = C„(0) . (6) 

It is easily seen that, via this change of variables, (5), the equations of motion 
([1]) satisfied by the quantities CniT) entail the following (autonomous) equations of 
motion (in the real time t) for the particle coordinates Zn{t): 

. , ■ gn+2 , 5n+l /-X 

Zn+liLlZn = \ . (7) 

Zn Zn-\-l Zji Zn+2 

Here and hereafter superimposed dots indicate differentiations with respect to the 
time t. 

So, this model ([7]) describes the "physical evolution" which we study. Note 
that its equations of motion, ([7]), are of Aristotelian, rather than Newtonian, type: 
the "velocities" i„, rather than the "accelerations" z„, of the moving particles are 
determined by the "forces" . In Appendix D we discuss the connection of this model 
with more classical many-body problems, characterized by Newtonian equations of 
motion. 

Let us immediately emphasize two important qualitative aspects of the dynamics 
of our physical model ([7]). The "one-body force" represented by the second term, 
iujZn, in the left-hand side of the equations of motion ([7]) becomes dominant with 
respect to the "two-body forces" appearing in the right-hand side in determining the 
dynamics whenever the {complex) coordinate z„ of the n-th particle becomes large (in 
modulus). Hence when |2:„(i)| is very large, the solution z„(t) of ([7]) is approximated 
by the solution of i„ + i z„ « implying that z„ [t) is characterized by the behavior 
Zn{t) ~ c exp (— i /;), therefore the trajectory of the n-th particle tends to rotate 
(clockwise, with period 2 T) on a (large) circle. This effect causes all motions of our 
physical model, d?]), to be confined. Secondly, it is clear that the two-body forces (see 
the right-hand side of ([7])) cause a singularity whenever there is a collision of two (or 
all three) of the particles as they move in the complex z-plane, and become dominant 
whenever two or three particles get very close to each other, namely in the case of 
near collisions. But if the three particles move aperiodically in a confined region (near 
the origin) of the complex z-plane, a lot of near collisions shall indeed occur. And 
since the outcome of a near collision is likely to be quite different depending on which 
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side two particles scatter past each other - and this, especially in the case of very 
close near collisions, depends sensitively on the initial data of the trajectory under 
consideration - we see here a mechanism complicating the motion, indeed causing 
some kind of chaos associated with a sensitive dependence of the motion on its initial 
data. This suggests that our model in spite of its simplicity, is likely to be rich 
enough to cause an interesting dynamical evolution. We will see that this is indeed 
the case. But before proceeding with this investigation let us interject two remarks 
(somewhat related to each other). 

Remark 1 This system ^ is still invariant under translations of the independent 
variable t (indeed, it is again autonomous j but, in contrast to (QP, it is no longer 
invariant under translations of the dependent variables Zn(t) nor under a simple 
rescaling of the independent variable t and of the dependent variables Zn{t). 

Remark 2 The general solution of the equations of motion ^ has the structure 

Zn{t) = ZcAlit) + Zn{t), (8 a) 

where the three functions i„(t) satisfy themselves the same equations of motion ^ as 
well as the additional restriction 

zi{t) + Z2{t) + Z3{t) = (86) 

which is clearly compatible with these equations of motion, and correspondingly zcM{t) 
is the center of mass of the system 

zijt) + z^jt) + zsjt) 
zcKiit) = ^ , (9 a) 

and it evolves according to the simple formula 

zcAiit) = zcm{0) exp{—iujt) — Z exp{—iujt) . (96) 

In Section [3] (and Appendix A) we determine the equilibrium configurations of 
our physical model, namely the values z!^'^'' of the three particle coordinates z„ such 
that 

0„ = z^-', i„ = (10) 

satisfy the equations of motion and we ascertain the behavior of our system in 
the neighborhood of these configurations. In the second part of Section [H and then 
almost always in the rest of this paper (and throughout the rest of this section) we 
restrict for simplicity consideration to the semisymmetrical case, see A main 

finding in Section [3] (and Appendix A) is that in the semisymmetrical case our model 
([7]) possesses generally two equilibrium configurations 2^°'''. We moreover determine 
the three exponents 7*-'"-' characterizing the small oscillations of our system in the 
neighborhood of each of these two configurations, defined according to the standard 
formulas (see Section IH]) 

z„(t) = z^°^'+eu;„(t), (11a) 

wi"'\t) = exp(-z7(")c^i)^'i"\ 

where of course e is an infinitesimally small parameter and the quantities wl™^ are 
time-independent. We find that the first two of these three exponents take in both 
cases the simple values 

7^=1, 7^') =2; (12) 
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the first of these corresponds of course to the center-of-mass motion, see (|9&|. As for 
the third exponent 7^^-', we find for one equilibrium configuration 

f + ^ M 



and for the other 



35 1 - M 

Here we have introduced the constant /i, 

whose value, as we shall see, plays an important role in determining the dynamical 
evolution of our model: in particular, this evolution does largely depend on whether 
or not /z is a real rational number, and if it is rational, 

(15) 

with p and q coprime integers (and q positive, q > 0), on whether the two natural 
numbers \p\ and q are large or small. A hint of this is already apparent from the 
results we just reported: while the solutions Wn\t) and Wn\t), see (11), are both 
periodic with period 2 T (see ()5c|) ; in fact Wn (t) is periodic with period T) , the 
solution Wn \t), see (11), is periodic with the period T, 
IT 

which is clearly congruent to Tonly if ^ is rational, see (jl5p and (13) - implying then 
that the small oscillations around the equilibrium configurations are always completely 
periodic with a period which is a finite integer multiple of T . 

In Section [3] we also introduce the special class of [exact and completely explicit) 
"similarity" solutions of our equations of motion, ([7]), and analyze their stability, 
namely the solutions of our system in the immediate neighborhood of these similarity 
solutions. 



2.3. Conserved quantities 



It is important to note at this point that the auxiliary model ([T]) possesses conserved 
quantities, which will be used in Section[5]to obtain its general solution by quadratures. 
Firstly, due to the translational invariance it is obvious that the quantity 



n=l 



(17) 



does not depend on r. In addition, the analysis of Section [5] shows that in the 
semisymmetrical case there exists an extra conserved quantity given by 



A'=(2C3-Ci-C2) 



1 - 



(Cl-C2)^ + (C2-C3)^ + (C3-Cl)^ 

2fi (2C3-C1-C2)' 



(18) 



Here the constant /i is defined in terms of the coupling constants g and /, see ([2]), by 
(|14p. We already mentioned that the value of this parameter (in particular, whether 
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or not is a rational number) plays an important role in determining the dynamical 
evolution of our model. A hint of this is now provided by the appearance of this number 
yU as an exponent in the right-hand side of ([T5)) . since this exponent characterizes the 
multivaluedness of the dependence of the constant K on the coordinates C„. 



3. Equilibrium configurations, small oscillations and similarity solutions 
of the physical model 

In this section we discuss, firstly, the equilibrium configurations of our physical model, 
([7]), and its behavior near equilibrium, and secondly, a special, explicit "similarity" 
solution of our model and its stability. 

The equilibrium configurations of our physical model (O, 

z„(t) = z^-', i„(0-0, (19) 

(see (ITOl) ) are clearly characterized by the algebraic equations 



^(oq) _ 



9n+l , 9n+2 



.(oq) _ (cq) 



(20) 



These algebraic equations entail 

^(cq) _^ _^(cq) ^ ^(cq) ^ Q _ ^^^^ 

It is now convenient to set 

z<,-' = (2 ^ Lo)-^'^ a„, (22) 
so that the equilibrium equations ([^0)1 read as follows: 

O^n _ grt+l _^ gn+2 ^23) 

2 a„ - a„+2 OLn - an+1 
These algebraic equations can be conveniently (see below) rewritten as follows: 

OLn = /3„+l {an - a„+2) + Pn^2 (^n " ttn+l) , (24a) 

via the position 

= -2 ■ (246) 

(Q!„_i - a„+i) 

We now note that, in order that the three equations (|24al) (which are linear in the three 
unknowns «„, although only apparently so, see (|246[) ) have a nonvanishing solution, 
the quantities /3„ must cause the following determinant to vanish: 

/32 + /33 - 1 -/33 -/32 
-/33 /33 + /3i - 1 -/3i 
-/32 -/?! /3i + /32 - 1 

To analyze the small oscillations of our system ([7]) around its equilibrium 
configurations we now set 

Zn{t) = z'^-'> +ew„{t), (26a) 

(see (|llap ) and we then get (linearizing by treating e as an infinitcsimally small 
parameter) 

Wn+iujWn+iUJ /3„+i (w„ - Wn+2) + Pn+2 (Wn " W^+l) = . (26&) 



. (25) 
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Therefore the three exponents 7^™) characterizing the small oscillations around 
equilibrium via the formula 

u;!") it) = exp(-z 7^'") lu t) w^™) , (27) 

providing three independent solutions of the system of linear ODEs (j266|) . are the 
three eigenvalues of the symmetrical matrix 

( P2+k + l -133 -h \ 

B = -/33 Pz + Pi + l -Pi , (28) 

V -(32 -Pi Pi +(32 + 1 J 

and the three 3- vectors w'-™-' = (^v^\ V2"^\ v^^^ are the corresponding eigenvectors, 

^i?„,x;(™)=7M„M . (29) 
£=1 

Hence the three exponents 7^™^ are the three roots of the "secular equation" 

P2+P3 + I-I -P3 -P2 

-P3 P3+P1 + I-"/ -Pi = . (30) 

-P2 -Pi /3i + /32 + 1 - 7 

Clearly these three roots are given by the following formulas: 

7''^-!, 7^'^ -2, 7(3) ^2(/?i+/32+/33) . (31) 

Indeed the determinant pop vanishes for 7 = 7*^^) = 1 (when each line sums to zero) 
and for 7 = 7^^^ = 2 (see (US])), and the third solution, 

7(3) =2(/3i+/32+/33), (32) 
is then implied by the trace condition 

trace [B] = 3 + 2(/3i + P2 + P3) = 7^^^ + 7^^^ + 7^^^ • (33) 
The first of these 3 solutions, 7^^^ — 1, corresponds to the center of mass motion (it 
clearly entails v^^ = v'^^\ see (gT]) and 

In the semisymmetrical case ^ the equations (or equivalently (24)) 

characterizing, via (j22p . the equilibrium configurations can be solved explicitly (see 
Appendix A) . One finds that there are two distinct equilibrium configurations (in fact 
four, if one takes account of the trivial possibility to exchange the roles of the two 
"equal" particles with labels 1 and 2), the first of which reads simply 

z<=<" = 0, zf^' = -z^^"' = (z"=^') ' = ^-J^^ (34) 

while the second has a slightly more complicated expression (see Appendix A). Note 
however that, in both cases, there holds the relation 

(zf" - z<-')' + (4°" - 4°")' + (4°" - ^r')' = ^ ■ (35) 

Moreover, in both cases the corresponding values for the eigenvalue 7^'^'', see are 
easily evaluated. The first solution yields (see (|I3ap ) 

(3) / + 8 g 1 9 ,or ^ 

7^ ' = r , o = - = -> 36a) 
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where, for future reference, we expressed 7''^-' not only in terms of the parameter fi, 
see ([H]), but as weh in terms of its rational expression (whenever applicable), 
while the second solution likewise yields 

T ' = — = -j = ■ (36&) 

Note that this implies that in the "integrable" case ^ both these formulas, l\'S6a\i and 
(|366p . yield 7'^' = 3; but it is easily seen that in this case only the first equilibrium 
configuration ([M]) actually exists. So in the "integrable" case the oscillations around 
the (only) equilibrium configuration (|34p are the linear superposition of three periodic 
motions (see ((27)) ) with respective periods 2T, T and ^ (see ([5c|)). Also in the 
"two-body" case (4) the second equilibrium configuration does not exist, while the 
first formula, (|36ap . yields 7^^-' = 1, so in this case the small oscillations around 
the (only) equilibrium configuration (|34p are the linear superposition of two periodic 
motions, with periods 2T and T (see ([5c|); consistently with the explicit solution, 
easily obtainable from (|¥c|) via (5)). 

As can be easily verified, the equilibrium configurations (|19p with ((20)) are merely 
the special case corresponding to ^c'Af(O) = 0, c = of the following two-parameter 
family of (exact) "similarity" solutions of our equations of motion ((7)): 

Zn{t) = zcM (t) + Zn{t; c), (37a) 

f„(i;c) = z^""' [1-Hc exp(-2^tJ^)]^/^ (37&) 

with the center of mass coordinate zcM{t) evolving according to ((9&)). The two 
arbitrary (complex) constants featured by this solution are of course zcniO) = Z 
(see ((Ml) and c, while the constants z'^'^^'s are defined as in the preceding section, see 

Clearly these (exact) solutions correspond, via the trick (5), the relation ((^^ 
(which is clearly consistent with ((^5)) and ((^0)) ) and the simple relation 

n = (38) 

to the two-parameter family 

C„(t) = Z + a„ (T-rb)l/^ (39) 

of (exact) solutions of (|T]). 

Let us now discuss the stability of this solution, ((37&p . To this end we set 

Zn(t) ^ Zn(t;c) + ewn(t), (40a) 

and we insert this ansatz in our equations of motion ([7]) , linearizing them by treating 
£ as an infinitesimally small parameter. We thus get 

: , . ~ , iu [Pn+l [Wn- Wn+2) + Pn+2 (Wn- Wn+l)] 

Wn+li^Wn^ — — — =0, (40&) 

1 -I- c exp (—2iLut) 



having used the definition ((24 &|) . Clearly the solution of this system of ODEs reads 
Wn(t) = exp(-iujt) Xn('&) , (41a) 

with 

I + c exp (—2iLot) 



1? = i9(t) ^t-(2iLu) ^ log 



1 



(416) 
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and the functions Xn (^) solutions of the autonomous linear system of first-order ODEs 

X'n+^^ iXn " Xn+2) + Pn+2 {Xn - Xn+l)] = 0, (41 c) 

where the primes denote of course differentiation with respect to i^. Hence (see (126 bp ) 
the three independent solutions of this linear system are 

Xl™)(^9)=exp(zc^^) u;(r)(z9), (41 rf) 

with the functions wi™'' defined by ((27)) (of course with t replaced by yielding 
via (P7|) and (|41ap with (j41bl) the following two equivalent expressions for the three 
independent solutions of the linear system (|40&[) : 



1 + c cxp (—2 iiut) 



1 + c 

(7<"'-l)/2 



(7'"'-l)/2 

exp{ij^"'^ujt)dl^''\ (42a) 



exp {2iLut) + c 



1 + c 



exp(it^i)Si'") . (426) 



Here the three exponents 7'™^ are defined as above, see (|3T|) . and likewise the 
"eigenvectors" wi™'' coincide with those defined above up to {arbitrary) normalization 
constants c'"*-', 

Note the equivalence of the two expressions (|42ap and (j426p (the motivation for writing 
these two versions of the same formula will be immediately clear) . 

For m = 1,2,3 these solutions, see (j42&[) . are periodic functions of the (real) 
time t with period 2r if |c| > 1. If instead |c| < 1, the solutions (see (|42ap ) with 
m = 1 respectively m = 2 are periodic with periods 2T respectively T (see ([3T|) ). 
The solution with m = 3 is periodic if 7'-'^) is real, but with the period ■pjjy which 

is not congruent to T if 7^'^^ is irrational; it grows exponentially with increasing time 
if Im [7'^'] < 0, implying instability of the solution (37) in this case, and it instead 
decays exponentially if Im [7^'^-'] > 0, implying a limit cycle behavior in configuration 
space, namely asymptotic approach to a solution completely periodic with period T or 
2 T depending whether the center of mass of the system is fixed at the origin or itself 
moving with period 2 T; but note that in this paper we restrict our attention to the 
case with real coupling constants. 



4. Analytic structure of the solutions of the auxiliary model 

In this section we discuss the properties of analyticity as functions of the complex 
variable r of the solutions Cn(''') of the auxiliary model ([l} (with arbitrary values of 
the 3 coupling constants gn, i- e. not restricted by the semisymmetrical condition 
1^: except when this is exphcitly specified, see below). In particular we show first of 
all that, for appropriate initial data characterized by sufficiently large values of the 
moduli of all three interparticle distances, namely by the condition (see ^) that the 
quantity 

Cn.in= min |Cn(0)-C™(0)| (43) 

n,m— n^m 
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be adequately large, the solutions ('«(''") a-re holomorphic in a disk Dq of (arbitrarily 
large) radius do centered at the origin, r = 0, of the complex r-plane (of course the 
"adequately large" value of the quantity Cmin depends on do, and on the magnitude 
of the three coupling constants gn', see (|5ip below). We moreover discuss via a local 
analysis a la Painleve the nature of the singularities of the solutions Cn{T) of the 
auxiliary model ([T]) as functions of the complex variable t and we thereby justify the 
assertions made in this respect in Section 2. 
To prove the first point, set 

^n(T) =Cn(T)-Cn(0), (44a) 

SO that these quantities cr„(T) vanish initially, 

a„(0) = 0, (446) 
and, as a consequence of ([1]), satisfy the equations of motion 



~ UO) ~ Cn+2(0) + a„(r) - a„+2(r) 

Cn(0)-C«+i(0) + 'T„(r)-a„+i(r) ■ ^ 

A standard theorem (see, for instance, [20j) guarantees then that these quantities 
crnir) - hence as well the functions Cn(''")i see ()44aP - are holomorphic in r (at least) 
in a disk Dq centered at the origin t = in the complex T-plane, the radius do of 
which is bounded below by the inequality 

(this formula coincides with the last equation of Section 13.21 of ^20], with the 
assignments m = 3 and a = oo, the first of which is justified by the fact that the 
system (|44cP features 3 coupled equations, the second of which is justified by the 
autonomous character of the equations of motion (|44c|) l. The two positive quantities 
b and M in the right-hand side of this inequality are defined as follows. The quantity 
b is required to guarantee that the right-hand sides of the equations of motion (j44cp 
be holomorphic (as functions of the dependent variables a„) provided these quantities 
satisfy the three inequalities 

k„|<&; (46) 

clearly in our case a sufficient condition to guarantee this is provided by the single 
restriction 

6<^, (47) 



with ^min defined by The second quantity in the right-hand side of ([iS]) . 

M = M{b), is the upper bound of the right-hand sides of (|44cP when the quantities 
(T„ satisfy the inequality (|46l) ; but of course the inequality ((45)) holds a fortiori if we 
overestimate M, as we shall presently do. Indeed clearly the equations of motion ()44cP 
with (gel) and (gZl entail 

AG , , 

M<-— ^, (48) 



with 



G= max |5„| . (49) 

n— 1,2,3 
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Insertion of 



in (gS]) yields 

b (Cmi„-2 6) 



do > 



16 G 



(50) 



hence, setting b — (to maximize the right-hand side; note the consistency of this 
assignment with (|T7)) ). 

Sn 



do > 



^2 

^min 



(51) 



128 G" 

confirming the assertion made above (that do can be made arbitrarily large by choosing 
Cmin adequately large). 

Next, let us show, via a local analysis a la Painleve, that the singularities as 
functions of the complex variable r of the general solutions Cn{T) of our auxiliary 
model ((!]) associated with a coincidence of two of the three components are square- 
root branch points (recall that a singularity at finite r of a solution Cn(''") of the 
evolution equations ([T|) may only occur when the right-hand side of these equations 
diverges). Such a singularity occurs for those values Tb of the independent variable t 
such that two of the three functions C„ coincide, for instance 

Clin) ^ C2in) ^ Csin) ■ (52) 

The square-root character of these branch points is evident from the following 
ansatz characterizing the behavior of the solutions of ([T|) in the neighborhood of these 
singularities: 



Ur) = Cb- i-1)' « (r - Tfc)^/' + Vs{t- n) + 



fc=3 



n) 



k/2 



n) + ^OL. 

k=3 



(fe) 
3 



(t - n) 



k/2 



with 



93, V3 



91 +92 



5 9s+i 



r ^ ' - P. (t t V s = l,2mod(2), 

Cb - C3& 6 (Cb - C3fc) 

and the constants (Mf' determinable (in principle) recursively (for fc = 
inserting this ansatz in ([1]), so that, to begin with 

2 a (52 - 9\) 



„(3) 



3 (Cb - C3)' 

- (-1)^ 



36 (Cb 
= 1,2 mod(2). 



C3)^ 



3 (ffs - 7 5s+i) + 



(ffi -92) 

93 



= 1,2 
(53 a) 

(53&) 
3,4,...) by 

(53c) 
(53rf) 



and so on. The diligent reader will verify the consistency of this procedure, for any 
assignment of the three constants rfc, Cb, and Csfc, which remain undetermined except 
for the obvious restrictions Tf, ^ 0, Cb 7^ 0, C3b 7^ Cb- The fact that (53) contains 
three arbitrary (complex) constants - the maximal number of integration constants 
compatible with the system of three first-order ODEs ([T]) - shows that this ansatz is 
indeed adequate to represent locally, in the neighborhood of its singularities occurring 
at T = Tb, the general solution of (fT|). 
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An analogous analysis of the behavior of the solutions of the system ([T]) near the 
values of the independent variable r where a triple coincidence of all three functions Cn 
occurs (corresponding to the excluded assignment (^^j, = in the above ansatz (53)), 
indicates, somewhat surprisingly, that such a triple coincidence, 

Clin) ^ C2{n) ^ Csin) ^ Z (54) 

might also occur for the general solution of the system ^ . This conclusion is reached 
via a local analysis analogous to that performed above, and is then confirmed (for 
the semisymmetrical case, see ([2])) by the exact treatment of Section [5l Indeed the 
natural extension of the above ansatz (53) characterizing the behavior of the solutions 
of ll]) in the neighborhood of such singularities, corresponding to a triple coincidence, 
see (|54l) . of the three functions Cn{T), reads as follows: 

C„(t) =Z + rj„iT- n) (1-^'/' + a„ (r - n) ^'^ + o (|r - n\ ^'^) , (55a) 

provided 

Re(7)<0. (55 &) 

Here the three constants a„ are determined, as can be easily verified, just by the three 
nonlinear algebraic equations ((23|) that were found in the preceding section while 
investigating the equilibrium configurations of our physical system ([7]) , while the three 
constants r/„, as well as the exponent 7, are required to satisfy the algebraic equations 

(7 - ^)'nn ^ 9n+l{rin ^ V71+2) ^ .9«+2(??n - V71+1) ^^g^ 

2 (a„ - a„+2)^ {an-an+i f 

These algebraic equations, and ([55)) . can be conveniently rewritten as follows: 

aji = /3n+i {an - an+2) + /3n+2 (^n - ^n+l) , (57a) 

(7 - l)r]n = Pn+1 (rin - f?n+2) + Pn+2 {Vn - Vn+l) , (57&) 



via the introduction of the quantities /?„ , see (j24&|l . Note that in this manner these two 
sets of equations, (|57a|) and (j576|l . have a quite similar look, which should however 
not mislead the reader to underestimate their basic difference: the three equations 
(|57ap are merely a convenient way to rewrite, via the definition (j246p . the three 
nonlinear equations (|23p . which determine (albeit not uniquely, see Appendix A) the 
three constants Q!„; while the equations (j57&P are three linear equations for the three 
quantities hence they can determine these three unknowns only up to a common 
multiplicative constant (provided they admit a nontrivial solution: see below). 

Of course these linear equations (j57b[) admit the trivial solution ?7„ = 0, and it is 
easily seen that there indeed is a special {exact) solution of the equations of motion ([T]) 
having this property, see ([39)) with the constants a„ determined by ([23)l and computed, 
for the semisymmetrical model, in Appendix A. This "similarity solution" (|39p of the 
system ^ has been discussed in the preceding section; but let us emphasize here 
that it only provides a iwo-parameter {Z and Tf,) class of solutions of the equations of 
motion ([1]), while the general solution of this system of three first-order ODEs must 
of course feature three arbitrary parameters. 

A general solution of the evolution equations ([1]) corresponds instead to the ansatz 
(|55al) if the linear equations (j57bp for the three coefficients rjn admit a nonvanishing 
solution, because in such mentioned above, a common scaling parameter 

for these three coefficients remains as an additional {third) free parameter (besides 
Z and Tf,). The condition for this to happen is the vanishing of the determinant of 
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the coefficients of these three hnear equations, (|576[) . namely again validity of the 
determinantal condition (|30p . a cubic equation for the unknown 7, which determines, 
as discussed in the preceding section, the three values (|3ip of this quantity. But 
the first two of these values, 7 = 7'^-' = 1 and 7 — 7'^^-' — 2 (see (I5T]) ). are not 
consistent with the requirement (j55bl) . The third solution, 7 — 7'^^ = 2(/3i + (32 + Ps) 
(see pip ) might instead be consistent with the requirement (j556p . and whenever this 
happens the ansatz (55) indicates that the general solution of the system of ODEs 
([T]) does feature a "triple coincidence" , see (f54| , and identifies the character of the 
corresponding branch point. 

In the semisymmetrical case ([2]) the equations characterizing the equilibrium 
configuration, ([23]) or equivalently (57), can be solved (see Appendix A). One finds 
that there are two distinct solutions of these nonlinear equations (j23p (in fact four, 
since each solution has a trivial counterpart obtained by exchanging the role of the 
two "equal" particles with labels 1 and 2). The first solution yields for 7 — 7^^^ the 
value p6ap . which is consistent with the condition (j55bp iff 

Re {^Ji) < 0, (58) 

and it yields for the branch point exponent, see (j55ap . the value 

while the second solution yields for 7'^^ the value (j36&p . which is consistent with the 
condition (j55&P iff 

Re(^)>l, (60) 

and it yields for the branch point exponent, see (I55ap . the value 

1-7 ^ M + 1 ^ p + q ^ 

2 2(m-1) 2ip-q)- > 

The last equality in ([59| and ([6T|) are of course only valid if /i is rational, /i = p/q. 

Note that these findings imply that the branch point associated with "triple 
coincidences" is not (only) of square-root type, being also characterized, see (j55ap . 
by the exponent the value of which depends on the value of the parameter /i, see 
(|59p and (|6ip : however this kind of branch point is not present if 

0<Re(/i)<l, (62) 

since in this case neither ([58)1 nor ([60)1 are satisfied. 

The results presented in this section are not entirely rigorous, since the local 
analysis of the singularities we performed above on the basis of appropriate ansdtze 
should be complemented by proofs that the relevant expansions converge. Moreover 
these analyses provide information on the nature of the branch points, but neither on 
their number nor their location. But these results are confirmed and complemented 
below (see Section [5]) by the analysis of the exact general solution of the equations 
of motion ([T]). Our motivation for having nevertheless presented here a discussion 
of the character of the singularities of the solutions of ([T]) via a local analysis a la 
Painleve is because an analogous treatment may be applicable to models which are 
not as explicitly solvable as that treated in this paper (see for instance [13] and [TT]). 
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In this section we obtain and discuss the general solutions of our models, ([T]) and ([7|). 
But since the general solution of the physical model ([7]) is easily obtained via the trick 
(5) from the general solution of the auxiliary problem ([1]), we focus to begin with on 
this model. 

A first constant of the motion is provided by the center-of-mass coordinate 

Ci + C2 + Cs 



z 



(63 a) 

(636) 
(63c) 



since the equations of motion ([T]) clearly entail 
Z' = 

hence 

Z{t) = Z{0) . 

And clearly the general solution of ([T]) reads 

Cn ir)^Z + Cn (t) , (64a) 

with the set of 3 functions C„ (r) providing themselves a solution of ([1]) , independent 
of the value of Z and satisfying the (compatible) constraint 

Ci (r) + C2 (r) + C3 (r) - . 
It is moreover clear that the equations of motion ^ entail 

CKi + C2 C2 + C3 Cs = .91 + .92 + .93, 
hence there also holds the relation 

Ci' + C2 + Cl = 2 (gi + 52 + 53) (r - To) . 
It is now convenient to set, as in the Appendix B of [13j . 



1/2 



p cos 



- (-1)^ -■ 

^ ' 3 



s = l,2. 



1/2 



p COS( 



Then, summing the squares of these three formulas and using the identities 

2 TT 2 TT 

cos(6') + cos(6' + — ) + cos(6' - — ) = 0, 

O O 
COS2(0) + COS^{e +—)+ COS^e -—)^-, 



one easily gets 
or equivalently 



(Cl-C2)' + (C2-C3)' + (C3-Cl)- 



hence, from (j656l 



P 



2 (<?i +92+ 53) (r - To) - 3 ^ 2 (51 + 32 + <73) (t - n) 



(646) 
(65a) 
(656) 

(66 a) 
(666) 

(67) 
(68) 

(69 a) 

(696) 

(69c) 
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3 

n=ro + — — — r, (69rf) 

2 (.91 + 32 + 93) 

which also entails 

P P = 91+ 92+ gz ■ (69e) 

Here we assume that the sum of the three coupling constants .g„ does not vanish, 
ffi + 52 + 53 7^ 0. The special case in which this sum does instead vanish is treated in 
Appendix C. The expression of the constant ri in terms of the initial data is of course 
(see dMcl) 

ri--.r^^ (70a) 

2(^1+32+53) ^ ' 



namely (see (|69cp ) 
or equivalently (see (|69dp i 



CK0) + Cl(0) + C|(0)-3Z^ 
2 (51 +52 +53) ' ^'"'^ 



(Cl-C2)' + (C2-C3)' + (C3-Cl)' . 

6 (51 + 52 + 53) 
There remains to compute 9 (r) , or rather 

u{t) = cos6'(t) . (71) 

Inserting the ansatz (66) in the equation of motion ([1]) with n = 3, one easily gets 

p2 (cose*)' (4 cos^ 9-1) = (4 31 + 452 + 53) cos6i - 4 (31 + 32 + 53) cos^ 

+ x/3 (31-32) sin . (72) 

From now on in this section - for simplicity, and because it is sufhcient for our 
purposes - we restrict attention to the semisymmetrical case ([2]), so that the last 
equation becomes simply, via ([71]) . 

p^u' (4«2_i) = (/ + 83) u-4 (/ + 2 3) . (73) 

The general case without the restriction ((2|) is treated in Appendix C. 

This ODE can be easily integrated via a quadrature (using ()69c|) ). and this leads 
to the following formula: 

1 



[u{t)] 



Ai-i 



^(r-ri), (74) 



where the parameter /i is defined by (jl4[) and K is an integration constant. Here 
we are of course assuming that f + 89 7^ (see (ITil) ): the case when this does not 
happen is treated in Appendix C. (Also recall that, as promised above, we shall treat 
in Appendix C the case in which the sum of the three coupling constants 3„ vanishes, 
namely when / + 2 3 = 0, which entails n = 0, see ([H])). As for the quantity K in 
(j74p . it is an (a priori arbitrary) integration constant. It is a matter of elementary 
algebra to express this constant in terms of the original dependent variables (via 
(|74p . (|69cP - ([TT]) and (66)), and one thereby obtains the relation 



K^l2{f + 2g)K (75) 

with K defined by (US]). This finding justifies the assertion that X is a constant 
of motion, see Section 2.3; and of course it determines the value to be assigned to 
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the constant K in the context of the initial-value problem. Likewise the value to be 
assigned, in the context of the initial-value problem, to the constant ti appearing in 
the right-hand side of ([71)1 is given by the formula 



(76) 



where (see ^ and ((666)) 1 



namely 



.(0, . ^<,m-<m-ao) 

2 { ICi (0) - <2(0)| ' + K, (0) - <3 (0)] = + (0) - <, (0)] = } 

Of course in these formulas the initial values Cn{0) of the coordinates Cn{T) of the 
auxiliary problem ([T]) can be replaced by the initial values Zn{0) of the physical problem 
O, see dll). 

Let us emphasize that we have now reduced, via (66) with (|69eP and ([711) . the 
solution of our problem ([T]) to the investigation of the function u{t) of the complex 
variable r, defined for t ^ as the solution of the (nondifferential) equation (|74p that 
evolves by continuity from u(0) at r = 0. 

To proceed with our analysis an additional change of variables is now convenient. 
We introduce the new (complex) independent variable ^ by setting 

( = (78) 

and the new (complex) dependent variable w = w {^) by setting 

w{0=^^^ [u{T)f . (79) 
Thereby the expression of the solution (66) of our original problem ^ reads 

1 /2 

s =1,2, (80a) 

1/2 

CAr)^Z-2(^l±^^ [^w{0]'^\ (80&) 

while the (nondifferential) equation that determines the dependence of w (S,) on the 
(complex) variable ^ reads 

[^(0-1]^-' [w(0]-''=i. (81) 

Note that this equation is independent of the initial data; it only features the constant 
12, which only depends on the coupling constants, see (|14p . 

We conclude that the solution of our physical problem ^ as the real time variable 
t evolves onwards from i = is essentially given, via (80) and (5), by the evolution 
of the solution w(^) of this (nondifferential) equation, (|8ip . as the complex variable ^ 
travels round and round on the circle S in the complex ^-plane defined by the equation 
(see dZH]) and ([5at ) 

^ = R exp(2iwi) R [exp (2 i w i) + 77] , (82a) 



Towards a Theory of Chaos Explained as Travel on Riemann Surfaces 



19 



namely on the circle with center ^ and radius \R\ . The parameters R and ^ (or 77) 
depend on the initial data according to the formulas (implied by (|78p . (5), (|75p. (fT5|) . 

R = Hf+M [1 _ ^] , (82&) 

2iu; [2z3(0)-zi(0)-Z2(0)]' 

^=RV, (82c) 

_ z c { [Ci (0) - C2 (0)] ' + [C2 (0) - C3 (0)] V [C3 (0) - Ci (0)] ' } 



3(/ + 2g) 
2/i [2C3(0)-Ci(0)-C2(0)] 



2 



[Ci (0) - C2 (0)] ' + [C2 (0) - C3 (0)] ' + [C3 (0) - Ci (0)] - 



(82d) 



(82e) 



Of course in these formulas the initial values C«(0) of the coordinates Cn(T) of the 
auxiliary problem ^ can be replaced by the initial values Zn{0) of the coordinates 
Zn{t) of the physical problem d?]), see ([6|). 

Let us emphasize that, as the complex variable ^ travels on the circle S - taking 
the time T to make each round, see (|82ap and ([5c|) - the dependent variable w {^) 
travels on the Riemann surface determined by its dependence on the complex variable 
^, as entailed by the equation ((8T|) that relates w (^) to its argument f - starting at 
t = from ^ — £_Q, 

^o = ^+R={ri + l) R, (83a) 

zujR {[c,{o)-C2 m ' + [C2 (0) - C3 (0)] ' + [C3 (0) - Ci (0)] ' I 

(see (82)) and correspondingly from w(^o) = wq, 

^ , 1 _ [Ci (0) - C2 (0)] ' + [C2 (0) - C3 (0)] ' + [C3 (0) - Ci (0)] ' .34. 

" « 2^^ [2(3(0) -Ci(0)-C2(0)]' 



(see (|82^ ). 

Let us therefore now discuss the structure of this Riemann surface, namely the 
analytic properties of the function w(^) defined by ((8T|) . There are two types of 
singularities, the "fixed" ones occurring at values of the independent variable ^, and 
correspondingly of the dependent variable w, that can be read directly from the 
structure of the equation ([81]) under investigation, and the "movable" ones (this name 
being given to underline their difference from the fixed ones) occurring at values of 
the independent and dependent variables, ^ and w, that cannot be directly read from 
the structure of the equation (j8ip under investigation (they "move" as the initial data 
are modified). 



5.1. Movable singularities 

To investigate their nature it is convenient to differentiate (I8ip . obtaining thereby 
(using again ((8T|) ) 

(85) 
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where the prime indicates of course differentiation with respect to ^. (Note that this 
ODE is impHed by the nondifferential equation ([5T|) . while its solution reproduces 
the nondifferential equation (I8f p up to multiplication of its right-hand side by an 
arbitrary constant). The position of the singularities, S,b, and the corresponding values 
of the dependent variable, Wb = w{£_b), are then characterized by the vanishing of the 
denominator in the right-hand side of this formula, yielding the relation 

Wb = (86) 

which, combined with ((8T|) (at ^ — ^b) is easily seen to yield 

^t^^l''^ ^r exp {2 71 ifik), fc=l,2,3, (87a) 



tb=Q -r exp 



2'Kpk 

q 



k^l,2,...,q, {87b) 



{^-l)-^[f^\ . (87c) 



In the last, (|87cl) . of these formulas it is understood that the principal determination 
is to be taken of the /i-th power appearing in the right-hand side. The first of these 
formulas, l|87aP , shows clearly that the number of these branch points is infinite if the 
parameter fi is irrational^ and that they then sit densely on the circle B in the complex 
^-plane centered at the origin and having radius r, see ()87c|) . Note that this entails 
that the generic point on the circle B is not a branch point (just as a generic real 
number is not rational); but every generic point on the circle B has some branch point 
(in fact, an infinity of branch points!) arbitrarily close to it (just as every generic real 
number has an infinity of rational numbers arbitrarily close to it). As for the second 
of this formulas, (j87b[) . it is instead appropriate to the case in which the parameter /i 
is rational, see (jlSp . in which case the branch points sit again on the circle B in the 
complex ^-plane, but there are only a finite number, q, of them. 

These singularities are all square root branch points, as implied by the following 
standard proof. Set, for ^ « ^h, 

w{0 = M + fl(e - ^b)^ + o (le - ^bf^''^^'>) , (88a) 
with the assumption (immediately verified, see below) that 
< Re {f5) < 1 . 



It is then immediately seen that the insertion of this ansatz in (|85p (is consistent and) 
yields 

1 2(1-^) _ ^ ^ A* 



Note that these results confirm the treatment of Section ID the square root branch 
points of w (f) identified here, see ([86|) . are easily seen to correspond, via (80), to the 
pair coincidence Ci(''fc) = C3{'''b) or C2{Tb) = C3(''fc); while there is an additional class 
of square-root branch points which only affect Ci(''") and C2(''')7 but neither C3(''") iior 
w {^) , and occur at 

w = 4^ (89a) 
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due to the vanishing of the second square-root term inside the curly bracket in the 
right-hand side of (|80al) . and correspond therefore to the coincidence CiC^b) = C2{Tb)- 
The corresponding values of ^ (as implied by (|89ap with (IHT|)) are 



(we use the plural to refer to these values because of the multivaluedness of the function 
in the right-hand side of this formula). 



5.2. Fixed singularities 

Next, let us consider the "fixed" singularities, which clearly can only occur at ^ = oo 
and at ^ = 0, with corresponding values for w. 

Let us investigate firstly the nature of the singularities at ^ = oo. Two behaviors 
of w(^) are then possible for ^ « oo, depending on the value of (the real part of) /x. 
The first is characterized by the ansatz 

u;(^) =a^ + o(|e|^^*''^) , Re(/3)<0, (90a) 

and its insertion in (|8ip yields 

/3 = --, = -exp(i7r^), (906) 
M 

which is consistent with (|90ap iff 

Re(/z) > . (90c) 
The second is characterized by the ansatz 

u;(0 = l + , Re(/3)<0, (91a) 

and its insertion in (|8ip yields 



(3=^—, 0^-1 = 1, (916) 
fi- 1 

which is consistent with (191 ap iff 

Re(/i) < 1 . (91c) 

We therefore conclude that there are three possibilities: if Re(/x) > 1, only the first 
ansatz, (90), is applicable, and it characterizes the nature of the branch point of w(^) 
at ^ = cxd; if Re(/i) < 0, only the second ansatz, (91), is applicable, and it characterizes 
the nature of the branch point of w(^) at ^ = oo; while if < Re(/^) < 1, both ansdtze, 
(90) and (91), are applicable, so both types of branch points occur at = oo. 

Next, let us investigate the nature of the singularity at ^ — 0. It is then easily 
seen, by an analogous treatment, that two behaviors are possible, as displayed by the 
following ansdtze: either 

w{0^aC^ + o{\^\^''^^^), Re(/3)>0, (92a) 
/? = --, = -exp(i7ryu), (926) 
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which is appHcable iff 

Re(/Li) < ; (92c) 

or 

w{0^l + a^^ + o(\^\^''^^^y Re(/3)>0, (93a) 

/3=^, a^-i = l, (93&) 
/i- 1 

which is appHcable iff 

Re{fi) > . (93c) 

This analysis shows that the function w{^) features a branch point at ^ = the nature 
of which is characterized by the relevant exponent f3, see ()92&P or (j936p . whichever is 
applicable (see (|92cp and (|93cl) '). But let us emphasize that there is no branch point 
at all at ^ = if neither one of the two inequalities (|92cp and (|93cl) holds, namely if 
< Re(^) < 1. 



5.3. Explicitly solvable cases 

Let us end this Section [5] by noting that the equation ([8T|) for certain rational values 
of fi reduces to such a low degree polynomial equation that it can be solved explicitly. 
In particular, the polynomial equation is of second degree ii fi = —1,1/2 or 2; it 
is of third degree if /i = —2, —1/2, 1/3,2/3,3/2, or 3; while it is of fourth degree if 
Ai = -3,-1/3,1/4, 1/2, 3/4, 4/3 or 4. 

The diligent reader might wish to use the corresponding explicit solutions formulas 
for these cases to verify the validity of the previous discussion. 



6. The physical model 

The solution (80) can also be written, via ([75]) . (|18p and (5), directly for the particle 
coordinates Zn{t), to read as follows: 

y.^t 2z3(0)-^l(0)-^2(0) , , „. ,.^,.1/2 

Zs[t) =Ze — — 77 exp (-2zwt) + 1 

6^ 

■(-[w{t)]"^ + {~r [l2^i~iw{t)]"^), s = l,2, (94a) 

l.\ 7 i'.ot 2Z3(0)-Zi(0)-Z2(0) , . „. , _,,l/2 r- ,,sil/2 

Z3{t) = Ze -^—^ [t] exp(-2iujt) + 1\ [w [t)\ , (94&) 

where the constant 77 is given in terms of the initial data by (|82(j|) and we set 

w{t) = w[^ (t)] , (95) 
so that this dependent variable is now the solution of the nondifferential equation (see 



(ED) 

[w {t) - 1] [iZ; (t)] = i? exp (2 i w t) + C = ^ [exp (2 i w t) + 77] , (96) 

where the constants R, ^ and r/ are defined in terms of the initial data, see (82) (and 
recall that the initial data Cra(O) can be replaced by the initial data Zn{0), see dH)). The 
dependent variable w{t) is of course the solution of this equation, ((96)) . identified by 
continuity, as the time t unfolds from t = 0, from the initial datum w (0) = wq assigned 
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at t = 0, see ([M)) : this specification is necessary, since generally the nondifferential 
equation (|8ip has more than a single solution, in fact possibly even an infinity of 
solutions. 

A discussion of the behaviour of this solution of the initial-value problem of our 
model d?]) with ([2]) clearly hinges on ascertaining how the solution w{t) of ((96)) evolves 
in time. This equation (|96p corresponds of course to the combination of (|8ip with 
(|82ap . Hence one must firstly elucidate the structure of the Riemann surface defined by 
the dependence ofw{^) on the complex variables ^ as determined by the nondifferential 
equations (|8ip . and then understand the consequences of a travel on this Riemann 
surface when the complex variable ^ evolves according to (j82ap . namely it travels 
round and round on the circle of center f and radius \R\ in the complex ^-plane. 

The first task is simple, its foundation being provided by the analysis provided 
in Sections [5 . 1 1 and \5l2\ 

The second task is much more demanding, inasmuch as it hinges on the detailed 
manner the sheets of the Riemann surface are connected via the cuts associated 
with the branch points discussed in Sections 15.11 and 15.21 The main results of this 
analysis have already been reported (without proofs) in [Tj; their derivation requires 
a sufficiently extended treatment to suggest a separate presentation [12]. To avoid 
unnecessary repetitions, also the detailed analysis of the Riemann surface associated 
to ([5T|) is postponed to [T^]. 

7. Outlook 

In this paper we report a deeper analysis of the model introduced in [1] explaining 
many results that were there reported without proof (such as the derivation of the 
general solution by quadratures) while adding some new material (such as a detailed 
analysis of the equilibrium configurations, small oscillations and similarity solutions 
of the model). The novelty of this approach in accounting for a new phenomenology 
associated to chaotic motion in dynamical systems lies in the fact that the solution 
is a multi-valued function of (complex) time, and a detailed analysis of its Riemann 
surface leads to very specific predictions in the simpler (/i rational) cases, while it also 
unveils a source of irregular behaviour in the more complicated (/x irrational) cases - 
unpredictable inasmuch as the determination of its evolution requires knowledge with 
arbitrarily large precision of the initial data. The full analysis of the dynamics of this 
model - including the geometry of the associated Riemann surface - is postponed to 
a future publication [T^] . 

The purpose of this series of papers together with other related projects [171 [19], is 
to go beyond the local analysis performed in the literature relating analytic properties 
of solutions in complex time with dynamical properties of the model (Painleve- 
Kowalewskaya and its non-meromorphic extensions) and to perform a full description 
of the global properties of the Riemann surface. This full description requires not just 
finding the type of branch points and their positions, but specifying how the different 
sheets of the Riemann surface are attached together at those branch points. 

Whenever possible, such an approach provides very detailed information on the 
dynamics that cannot be obtained by the more classical local analyses. 
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Appendix A 

In this appendix we solve, in the semisymmetrical case, see the nonlinear 
algebraic equations (j23[) that characterize the equilibrium configurations and we 
thereby compute the "eigenvalue" 7^"^-', namely we obtain its two expressions p6ap 
and ([M6|) . 

The equations to be solved read (see (|23l)) 

ai = 1 , (97a) 

ai — as ai — ol-2, 

a2 = , (97o) 

a2 — as ai — 012 

a3 = ^^ + ^^, (97c) 

a-i — a\ a^ — 0L2 

and they of course imply the relation 

ai + a2 + 013 = . (98) 
It is now convenient to set 

S = a.\+a.2, D = ai — a2, (99a) 

entailing 

S + D S — D 

ai = — - — , a2 ^ — - — , as ^ -S . (99&) 

From (the sum of) (|97ap and (|97&P we easily get 

S (9 - D^) ^2Ag S, (100) 

and from this we get two types of solutions. The first solution is characterized by 
S = 0, implying (see (|99b|) and ((g7a|) ) 

a3 = 0, ai = —a2 = a, ^ f + 2g, (101) 

entailing (via (I22p ) the solution ([34|) for the equilibrium configuration, as well as (via 
(j24&l) with ^) the expressions 

A=A = 7^, (102!,) 
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hence, via ([321), the first expression, (|36ap . for 7^^^. 
The second solution is characterized by 

9 52 - L'^ = 24 g . (103) 

We now subtract (|976p from (|97ap and we thereby easily get 

hence, via the preceding relation, 

^ 3 /, - ttll . (105) 

And via (j24&l) with ^ and (j996P this is easily seen to yield 

P,+f3^+P^^l±ll^ (106) 
6 g 

namely, via ([321), the second expression, (|36bp . of 7'^^^. Note moreover that, in both 
cases, one gets the relation 

(ai - 02) ^ + (a2 - as) ^ + (as - ai) ^ = 6 (/ + 2 g) , (107) 

as can be easily verified from poip as well as from ^ with pOSp . 



Appendix B 

In this Appendix we consider certain nongeneric (classes of) solutions of our physical 
problem ^ , characterized by special subclasses of initial data. 

If the initial data are such that 77, hence as well ^, vanish, rj — = Q - and this 
entails that the initial data satisfy the condition 

N(0) - ^2(0)] ' + [Z2(0) - Z3(0)] ' + N(0) - zi(0)] 2 = HI±1A^ (108) 

I OJ 

see (82) and ([HI); hence these initial data are not generic, depending only on 2 arbitrary 
complex parameters rather than on 3 such parameters (or, equivalently, only on 1 
rather than 2 such parameters besides the trivial constant Z that only affects the 
center-of-mass motion, see ([6])) - then the time evolution of the solution z„(i) of our 
physical problem see ([HJ, is clearly periodic with the period ^ rather than T. 
The consequence of this fact are sufficiently obvious not to require any additional 
elaboration. An example of this type is that characterized by the parameters 

Lo = f = 2n, g^n, ^ T = 1, M = - = ^, (109a) 

2 90 

and the initial data 

zl(0) = 0.2, Z2(0) = -0.79658 + 0.717791, 23(0) = 0.59658 - 0.71779i 

that imply that the center of mass is initially at the origin and therefore stays there 
for all time, Z = 0. These initial data are easily seen to satisfy the condition (|108p . 
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Appendix C 

In this Appendix we explain how to integrate the ODE ([72| in the general case when 
the three coupling constants gn are all different, namely when the restriction ([2]) 
identifying the semisymmetrical case does not apply, and we also provide the solution 
of the ODE ([72l) in the two special cases (belonging to the semisymmetrical class 
characterized by the restriction ([2])) the treatment of which had been omitted in 
Section [51 and as well in another special case not belonging to the semisymmetrical 
class. 



Solution of equation |7i?[ j in the general case 

In this subsection of Appendix C we indicate how the ODE ((72)) can be integrated 
in the general case when the three coupling constants (?„ are all different. It is then 
convenient to set 

y (r) = tan [6* (r)] , (110) 
so that this ODE reads 

V V {V^ - 3) 1 

(^2 + 1) (^AV^ + CV^ + AV + C-2) ^ (r-Ti) 

with 

^ ^ ^/3_(g]_-^2)_ 4gi+4g2 + g3 ^-^2) 

2 (ffi + .92 + .93) ' 2 (.91+32+53) ■ 

To integrate this ODE we set 

AV^ + CV^ + AV + C -2 = A {V -Vi) {V -V2) [V -V3) , (113) 

so that the three quantities Vn are the three roots of this polynomial of third degree 
in V. We then decompose this rational function of V in simple fractions. 



V (1/2 _ 3) 



5 



(V^ + l) {AV^ + CV'^ + AV + C-2) j-[V-Vj' 
where of course 

= i, V5 = ~i, (115) 
and the five quantities /ij are easily evaluated in terms of the 3 roots Vn'. 

5 

= VS- (-3 + {V,~Vk)-'. (116) 

k=i,k^j 

The integration of the ODE (|llip is now trivial (using (|114p ). and it yields (using 
(|115|) ) the final formula 

[V (r) - zr [V (r) + zr f[ [V{r) - K]^" ^ K {r - n) , (117) 
where K is the integration constant. 
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Solution of equation ^7S^ in two special subcases of the semisymmetrical case 

In this subsection of Appendix C we provide the solution of the ODE ([7^ in the 
two special subcases (of the semisymmetrical case) the treatment of which had been 
omitted in Section [5l and as well in another special case not belonging to the 
semisymmetrical class. If 

31+52+53-0, (118) 

p is constant (namely r-independent, p (r) = p(0), see (|69cp V Moreover, via the 
restriction ^ characterizing the semisymmetrical class, we get (see also 

f^~2g, (119a) 

Then ((74)) is replaced by 

"3/ (r-To)" 



, (r) exp [-2 (r)] = exp , (q) 
Let us also note that, if ([2|) were replaced by 



(119&) 



51 = -52 = 5, 93 = 0, (120a) 
which is also consistent with (|118p . then (fM)) with (|7ip would be replaced by 

^(r)+sm[2^(r)]^ ^^;,(;°--) . (120.) 

Returning to the semisymmetrical case characterized by validity of the restriction ([2]) 
we now consider the second case whose treatment had been omitted in Section \E\ 
namely 

/ = -85- (121) 
Note that in this case /i diverges, see Then ([74]) is replaced by 

u (r) exp [u^ (r)] = [K (t - rOr^/" . (122) 



Appendix D: relation with more standard (Newtonian) three-body 
problems 

In this Appendix we indicate the relation among the three-body problems treated 
in this paper, characterized by equations of motion of Aristotelian type ("the 
particle velocities are proportional to assigned external and intcrparticle forces"), with 
analogous many-body problems characterized by equations of motion of Newtonian 
type ( "the particle accelerations are proportional to assigned external and intcrparticle 
forces"). The results reviewed in this section arc of interest inasmuch as they relate 
the model treated in this paper to other, somewhat more physical and certainly more 
classical, many-body problems, including a prototypical three-body model introduced, 
and shown to be solvable by quadratures, by Carl Jacobi one and a half centuries ago 
|21| . and the one-dimensional Newtonian many-body problem with two-body forces 
proportional to the inverse cube of the intcrparticle distance introduced and solved 
over four decades ago (firstly in the quantal context [SJ [6] and then in the classical 
context 23 [26]), which contributed to the bloom in the investigation of integrable 
dynamical systems of the last few decades (see for instance [27, ) . 
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By differentiating the equations of motion ([T]) and using them again to ehminate 
the first derivatives in the right-hand sides one gets the following second-order 
equations of motion of Newtonian type: 



C" 

Sn 



{Cn — Cn+2)^ {Cn — Cn+l)^ 
{Cn — Cn+2) (Cn+2— Cri+l) (Cn " Cn+l) (Cn+1 — Cn+2) 



yn+2 



[gn - gn+i) 



(123) 



Likewise from the equations of motion ([7]) one gets 



2 5n+l '^gn+2 



gn+1 (gn - gn+2) I 

{Zn — Zn+2) {Zn+2 — ^n+l) (^n ~ Zn+l) (^n+l ^ Zn+2) 



gn+2 {gn — gn+i) 



(124) 

Of course the solutions of the first-order equations of motion, ((T|) respectively ([7]), 
satisfy as well the corresponding second-order equations of motion, p23p respectively 
(|124p . but they provide only a subset of the solutions of the latter. On the other hand 
it is again true that the solutions of the second-order equations of motion p23p and 
(|124p are related via the trick. 

In the integrable "equal-particle" case, see ([3]), these equations of motion simplify 
and correspond respectively to the Newtonian equations of motion yielded by the two 
standard iV-body Hamiltonians 

N 



g 



N 2 

n=l m,n=l;m=in ^ ^™ > 



respectively 



N 



2 ^2 



H [z, P)^2^ ^ 

n=l ^ 



N 



2 ( z^ ZjYi) 



(125) 



(126) 



m.n— 1; m^n 

with = 3, the complete integrability of which is by now a classical result (even in 
the A''-body case with > 3: see for instance [8]). 

In fact the more general three-body Hamiltonian models 



respectively 



H (C, 1l) 



3 

E 

n=l 



3 

E 

n=l 



2 



(Cn+l ~ Cri+2)^ 



pI + zl 



{Zn+l — Zn+2) 



(127) 



(128) 



featuring three different coupling constants that yield the equations of motion 



C" = - 



2 5n+i 2 gn+2 



(Cn ~ (Cn — Cn-l-l)^ 



respectively 



Z„ 4- W Zn 



2 ffra+l 2 gn+2 

{Zn ~ Zn+2)^ {Zn — Zn+l) 



(129) 



(130) 
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are also solvable by quadratures. For the equations of motion (|129p this discovery is 
due to Carl Jacobi [H] ; while the solutions of the equations of motion p30p can be 
easily obtained from those of the equations of motion (|129p via the trick (5) . For a 
detailed discussion of these solutions, and additional indications on key contributions 
to the study of this problem, the interested reader is referred to [13] and [8]. But we 
will perhaps also revisit this problem, because we believe that additional study of these 
models, (|129|) and (|130p . shall shed additional light on the mechanism responsible for 
the onset of a certain kind of deterministic chaos, as discussed above. 

Finally let us recall that our Aristotelian model ([7]), as well as the Newtonian 
models described in this Appendix, describing "particles" moving in the complex z- 
plane, can be easily reformulated as models describing particles moving in the real 
plane, with rotation-invariant (or at least covariant) real two-vector equations of 
motion (see for instance Chapter 4 of [8]). 
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